SPATIAL COAGULATION WITH BOUNDED COAGULATION RATE 
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Abstract. We prove that the spatial coagulation equation with bounded coagulation rate 
is well-posed for all times in a given class of kernels if the convection term of the under- 
lying particle dynamics has divergence bounded below by a positive constant. Multiple 
coagulations, fragmentation and scattering are also considered. 



1. Introduction 

This article is concerned with Smoluchowski's picture of a chemical reaction involving two 
kinds of dynamics. Particles involved in the reaction are described by their location x G M n 
and an inner structure, which we describe as an element y of an Abelian semi-group. It 
can be its mass, in M* or N, or a finer structure saying for instance which smaller masses 
compose that mass, as in |NorOO] . On the one hand, particles move as a result of the action 



of some convection flow and thermal diffusive effects. On the other hand, close particles may 
coagulate to form a new particle with a new characteristic. We shall mainly be interested 
here in binary collisions. Not to overload notations, let us describe a particle by its position 
iGl" and its mass y G 

The evolution of the system is described by a time dependent kernel (j,t(x ; dy) from M n to 
W + , with /if (x ; (a, /3)J representing the concentration of particles at location x with a mass 
in between a and /5. Let K be a measurable symmetric function on M*. x Kl. Given a signed 
measure \x on M^_, define a signed measure K(/j,,/j,) on setting 

K (n,fj>) = \ J {<W ~ s v~ s v'} K {y,y') v{dy)K d y')- 

Smoluchowski's picture of the microscopic dynamics leads heuristically to the formal evolu- 
tion equation on kernels 

(1.1) d t n t = L* x fi t + K(fi t ,fi t ), 

where L x — | ^(x, y) d Xi d Xj + b l (x, y) d x . is a convection-diffusion differential operator. This 
heuristic result was proved in an important special case by Norris in |Nor04] . and later 
in |HR07a| and |YRH09) by Hammond, Rezakhanlou and Yaghouti. Well-posedness for 



equation ( 11. II) has been investigated in a relatively small number of works. Amann's article 
|AmaOO] is the only work where the problem is tackled in this generality. Supposing all the 
coefficients in the dynamics bounded, he obtains the existence of a unique maximal solution 
depending continuously on the initial condition using highly non-trivial functional analytic 
tools. Under some smoothness and moment hypothesis on the initial condition, he is able 
to prove well-posedness for all times in space dimension 1 or if the diffusivity coefficients do 
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not depend on y, and b = 0. This work was partly extended in |AW05j . in the Brownian 
case where 6 = and 

(1.2) ^a ij (x,y)d Xi d Xj = d(y)A x , 

with d(y) bounded away from and oo, and well-posedness for all times was proved under 
some moment and size conditions on the initial condition. They allow linear growth for the 
coagulation rate. The works |Mis03| . |LM02b] . |LM02a] . |CDF10] . identify different settings 



in which one can prove the existence of a weak solution to equation ( II. ip (in different 
appropriate senses), defined for all times and for unbounded coagulation kernels controlled 
in different ways. They all investigate equation ( 11.11) without convection and with a diffusion 
of the form ( II. 2p . Hammond and Rezakhanlou are able to give in |HR07b] a well-posedness 
result for all times for a discrete mass space, a bounded decreasing diffusivity d(y) , and some 
moment conditions. (See also |RezlO| .) 

All these works investigate the differential problem (11.11) in different functional settings. 
Norris |Nor04] and Laurencot-Mischler |LM02bj consider a mild integral formulation of equa- 



tion (II. ip rather than its differential counterpart. They work with no convection term and 
a diffusivity of the form (11.21) . Relying crucially on Gaussian estimates, Norris is able to 
prove a well-posedness result for all times under some weak control on the coagulation rate 
and some moment condition on the initial condition. His method does not apply to position 
dependent diffusivity dynamics. Consult also [RezlO] for a recent existence and uniqueness 
result under some moment assumptions. 

We refer the reader to the above works for further and older references on the spatial 
coagulation equation (II. ip . Note yet that the methods developped in |AW05| and [LM02aj, 
|LM02bj . should be sufficiently robust to deal with more general dynamics than (II. 2p on 
bounded smooth domains as they essentially rely on some compactness property of the 
operator. 

Following Amann's view |Ama00| . we consider in this work a bounded coagulation kernel 
and a general convection-diffusion spatial dynamics. It is quite likely that under some 
reasonably general conditions on the convection-diffusion coefficients equation (II. ip is well- 
posed for all times in a good functional framework. Despite all the subtle works done, this 
remains yet an open problem. We identify in this work a previously unnoticed case where this 
problem can be solved. Roughly speaking, a condition on the divergence of the convection 
term prevents the particles in the underlying microscopic dynamics from concentrating too 
much; this enables us to obtain some uniform bound in the spatial coordinate. Although 
far from opening a road to solving the general problem, the simplicity of our framework and 
argument contrasts with the sophistications developped in the above references. 

We suppose that the coefficients a'- 7 and b % are measurable functions of (x, y) G M n x , of 
class C°° as functions of x G M n , and with all their derivatives bounded, uniformly in y G ; 
this allows linear growth of the drift. To give an integral form of ( II. ip . let us suppose here 
that one can write the matrix a as aa* for some bounded matrix-valued measurable function 
a on K n x of class C°° as a function of x. Denote by <p t {x,y) the well-defined random 
flow on R n x IR* of the stochastic differential equation 

dx t = a(x t , y t ) dB t + b(x t , y t )dt, 
{ - 3) dy t = 0. 
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For a kernel v on W 1 x M.* +J and t ^ 0, set 

(P t u) (x ■ dy) = E r e - Jo^X^v).!/)* I/ (0 t ( Xj y) ; 

where divfe is the divergence of the vector field b(-,y) on R n , for fixed y. The integral form 
of equation (II. ip is 

ft 



;i-4) fH = p tf jt + t P t _ s (K(fx s ,^ s ))ds. 

Jo 



We now introduce the functional setting in which we are going to solve this equation. Let 
is be a non-negative finite measure on such that 

(1.5) u -ku ^Ci/ 

for some positive constant C. Define M as the set of kernels v[x ; dy) from W 1 to of 
the form {f(x; y) ^(dy)}^ Rd , where f(x; y) is a bounded measurable function uniformly 
equicontinuous in x: 

Ve>0, 3J>0,V2/Gi;, Iz-x'l < <5 =^ \ f{x; y) - f(x'; y)\ < e. 

Setting 

= sup ^o-esssup \f(x; y)\ 



defines a norm on M. 

Theorem 1. Suppose the following two conditions hold. 

(i) There exists some positive constant a, /3 such that a Id ^ a(x, y) ^ Id, for all 
(x,y) G M n x E* . 

(ii) There exists a positive constant e such that the vector b(-,y) has a divergence bounded 
below by e, for all y G M^j_. 

Let M > be an upper bound for K. Then, for any fi G M , equation (jl.4j) has a unique 
maximal solution in N ' , continuous with respect to t, and defined for all times provided ||/i || 
is small enough; it is non-negative if /i is non-negative. 

Remarks. 

(1) Use Laplace transform to construct a positive continuous and LEB-integrable function 
g on such that the finite measure Vo(dy) := g(y) dy satisfies (II. 5p for some positive 
constant C. Suppose fio(x ; dy) = ho(x; y) dy, for a bounded measurable function ho 
with support in IR n x (0, m), for some m < oo. It follows from theorem [1] that \i t 
is well-defined for all times provided ||^o||oo is small enough, and each \i t {x\dy) is 
absolutely continuous with respect to Lebesgue measure. 

(2) It will be clear from the proof of theorem [T] that the conclusion remains unchanged 
if the coagulation rate K(y, y') is allowed to depend on x also, as an equicontinuous 
function with a modulus of continuity independent of y, y' G 

(3) Although all our statements are given here and below with y G M^_, all the results 
hold true with y taking values in a general Abelian semi-group describing the inner 
structure of particles, at the price of introducing the formalism used in |Nor0 0j. All 



results hold in particular for the discrete masses spatial coagulation equation. 
(4) Note that since we are working with a bounded coagulation kernel we do not require 
any moment assumption on the initial condition. 
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(5) A similar result holds when M. n is replaced by a bounded domain of M", with smooth 
boundary, and Dirichlet or Neumann conditions are imposed on the boundary Con- 
dition (ii) makes unnecessary the use of spectral properties of the diffusion/convection 
operators. 

(6) Our method of proof is similar in spirit to the method used in [AW05J, where the 
quadratic coagulation term is balanced by a linear termQ Our functional setting is 
more elementary, and the above result new. 

(7) The introduction of the map P t is better understood by considering first a pure 
convection dynamics, where a (and so a) is null. In that case, the equation [x t = L x fi t 
on kernels reads 

(1.6) (i t = -b % d Xi iH - (div6)/i t . 

Introducing the flow <p t of the ordinary differential equation x = b(x) on ~R n , equation 
(II. 6p has solution 

tH(x ; dy) = fioMx) ; dy) e fi(to">XM*),v)*> 

In the case of a general second order differential operator L x as above, the associated 
dynamics is no longer deterministic, but given by the stochastic differential equation 
(II. 3p . Equation (II. 4p is the integral form of (jl.ip obtained by varying the constant. 

2. Proof of theorem [TJ 

2.1. Well-posedness. The proof of the following lemma is straightforward and left to the 
reader. 

LEMMA 2. The space (A/", || • ||) is a Banach space. 

The proof of theorem [1] is an application of the elementary fixed point theorem for con- 
tractions of a Banach space. Fix T > 0. Given an A/"-valued path v = (v s )o^s^t, define for 
each t G [0, T] the kernel 

-i 



F(u) =P t uo+f P t - a (K(u a ,v a ))d8. 
Jo 



PROPOSITION 3. The map F is a locally Lipschitz map from C([Q,T],J\f) into itself. 

It follows classically from proposition |3] that one can associate to any fiQ G M a time T 
depending only on ||/io||, in an increasing way, such that F has a unique fixed point in 
C([0,T],A/") with initial value /xq. The following proposition proves then that T = oo. 

PROPOSITION 4. We have \\ji t \\ < for all t G [0, T], if ||/i || is small enough. 

We now proceed to proving the two propositions. 

PROOF OF PROPOSITION [3] - We first deal with the operator P t , for t > 0. 

LEMMA 5. a) The linear map P t sends J\f into itself and has norm no greater than e~ et . 
b) Fix v G M . The map t G M+ i-> Ptv G (ftf, \\ ■ ||) is continuous. 



1 The main ingredients of |AW05| are mass conservation of solutions and the exponential decay of a semigroup 
on functions with zero spatial mean. We do not need mass conservation but use crucially the exponential 
decay of Pt. 
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PROOF - a) For fi(x ; dy) = f(x; y)u (dy) in Af, 

(Ptfi)(x;dy) =E[e-J?^«^C»*^ d '/(0 t (x,y);y)] 

so we need to see that the above expectation defines a (y-) uniformly equi-continuous 
function of x; denote it by (P t Alvb f) (x, y). Rewrite equation (II. 3p as a stochastic 
differential equation driven by vector fields Vi 

dx t = Vi(x t , y) dB\ + b(x t ) dt, 

with Vi(x,y) being the i th column of the matrix a(x,y). Kusuoka and Stroock have 
proved in |KS87] that we have^ for each i G {1, . . . , n} 



f2.r 



sup IViP^fUx^) ^ ct'i sup \ f(x; y 



for each y G M.* + , where the constant c depends only on an upper bound for the uniform 
norm of the VJ's and all their derivatives. It can be taken to be independent of y G M^j_ 
from our assumptions. Also, from the uniform ellipticity of the matrices a(x,y), 
inequality (12. II) implied the (y-) uniform equicontinuity of P t dlvb f. This proves that 
P t sends M into itself; the statement on its norm is straightforward from hypothesis 
(ii) of theorem [TJ 

b) It suffices to prove the continuity at t — 0. For u(x, dy) = f(x ; y)h> (dy) in A/", we 
have 



(2.2) 



\P t u — u\\ ^ sup E 
s= sup E 
^ sup E 



e" fo^XM^y) Offafa y)-y)- f( X ; y) 

\f(M x ,y);y) - f( x ; v) 
\f(M x ,y)\y) - f( x ~i y) 



oo sup E 

x,y 



■ f*(divb)(<f> s (x,y) \y)ds _ ^ 



-t. 



for some positive constant c. By the uniform equi- continuity of f{-',y) one can asso- 
ciate to any a > a 5 > such that 

VyGi;, Vi.i'gR", \x'-x\ <<5=^ \ f(x,;y) ~ f(x';y)\ ^ a. 

Then, we have for any (i,n)6K"x M*_ 

E[\f(4> t ( x ,y);y) - f( x ;y)\] < aF(\cf> t ( x ,y) ~ x \ < 5)+2\\f\\ 00 F(\<l> t (x,y)-x\ > 5) 

E[|<Mx,2/)-x| 2 ] 



^ a + 2 



( 9 



They prove their results for a and b and all their derivatives bounded; their proof of the following result 
extends in a straightforward way to our elliptic setting. 

3 Due to the uniform ellipticity of the matrices a(x,y), inequality (|2.ip is equivalent to the existence of a 

constant c', independent of i 6 {1, n}, such that sup |<9 a;i P t dlvb /| ^ c't~ x / 2 sup |/(a;,y)|. The y-uniform 

igi" igi" 
equi-continuity of p^ lvb f follows from the boundedness of / as t > is fixed. 
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As the isometry property of stochastic integration with respect to Brownian motion 
and Jensen's inequality give 



E[\(P t (x,y)-x\ 2 ] ^2E (/ a(<i> a (x,y),y)dB. 



< 2E 



\a((f) s (x,y),y)\ ds 



'-Jo 



+ 2E 
2tE 



1 X2 

b((p s (x,y),y)ds 



o 



\b((j) s (x,y),y)\ ds 



'-jo 



^ Ct 



for some positive constant C depending only on \\a\loo and 



sup 



E 



\f(M x >y) ;y) - f( x ; v)\ 



< a 



2C\\fl 



5 2 



it follows that 



t. 



Together with (12. 2p . this inequality implies that \\PtV — v\\ — > as t goes to 0. □ 

It follows from lemma that proposition [3] will be proved if we can prove that the kernel 
K(fi, /i) belongs M if \x does (c), and is a locally Lipschitz function of \x (d). Recall we 
write M for an upper bound of K. 

c) Given n(x, dy) = f(x ; y) vo(dy) in Af, set, for any x G M n , 

1 



AT + (/U, fi)(x ; dz) 



5y +y ,(dz) K(y,y') f(x;y') vo{dy') f(x; y)v {dy), 



and 



K {p,,n)(x;dz) = yj K(y,z)f(x;y)u (dy)Jf(x;z)v (dz), 

so that K(/j,, //) = — K~(/J,,fi). As z^o satisfies hypothesis (jl.5p we havd^ 

< Af||/||> ^o ^ CM\\f\\lu ; 

we can thus write 

K + (fx, fj,)(x ; dy) = k + (x, y) u (dy) 

for some bounded measurable function k + (x, y). Taking any bounded function g, we see 
by dominated convergence that the integral 

g(y) K + (fi,[i)(x;dy) 

is a uniformly continuous function of x, with a modulus of continuity independent of g. It 
follows that one can construct a version of k + (x, y) which is (y-) uniformly equicontinuous 
in x. So K + sends M into itself; it is straightforward to see that K~ does the same. 

d) Denote by | • \jy the total variation norm on finite signed measures on R+. The 
elementary inequality 



(K(fi,fi)-K(fi',fi'))(x; -)\ TV ^ m(\h(x]-)\tv+\v'(x]-)\tv) \ fi(x ; ■)-(!' (x ; ■) 
shows that K is locally Lipschitz in (A/", || ■ ||). 



TV 



x G 



We write \p\ for the absolute value of a signed measure p. 
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PROOF OF PROPOSITIONS! - Set m = u (M* + ). As equation (fOj) gives for all t G [0,T] 

H/Xtll «S e- e *||// || +M(^ + m) jf e- £( '- s) ||^frfs, 
the real- valued function \\pt\\ is bounded above by the solution of the equation 



z t ~- 
Its solution 



ez t + M^ + rnjzl z Q = ||// ||, < t < T. 



n r/r \ e-Af(£+m) 

M(% + m) H v I, o 



w " et 



is defined for all times and less than or equal to Zq for \\uo\\ < — — t> 

2.2. Positivity. We do not suppose in this section that /j,q G J\f is small, so we work on a 
time interval [0, T] where equation (jl.4p is known to have a unique solution in J\f continuous 
with respect to t. 

Suppose fj,Q ^ 0. We prove the positivity of p t by the usual method; see for instance 
theorem 10 in |AW05] . We recall it here for the reader's convenience. Let a be a positive 
constant such that 

(2.3) / /^(x;%)^|HK(R;)^ 

Jr* + m 

for all x G and s G [0, T]. Define a map G from C([0,T],J\f) into itself setting 

G((u s ) 0<s<T ) t = P t p + / P t - s (K(fi s , u s ) + av s ) ds- a / P t - S u s ds. 

Jo Jo 

Given two measures p and p' on set 

K(p,p')(dz) = J {5 y+y ,(dz) - S y (dz) - 5 y ,(dz)}K(y,y') p(dy)p'(dy'). 
As the following inequality between kernels holds pointwise in x G lR n 

K(fj, s , v s )(dz) ^ -( / K(y,z) p s (dy))u s (dz) + au s (dz) 

>{a- / K(y,z)fi s (dy))u s (dz) > 0, 

the first integral in the definition of G is non-negative. The path (/is)os;sigT solves by con- 
struction the equation on kernels 

u t + a [ Pt-sV s ds = G((v s ) 0<s ^ T ) t >0, te[0,T]. 
Jo 

Taking T smaller if necessary, one can solve this equation by a fixed point method; see e.g. 
theorem 10 of |AW05j . As each iteration step of the method preserves non-negative kernels, 
its unique fixed point (p* s )o^s<:T is non-negative if /xq is non-negative. 
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REMARK ON THE PROOF OF THEOREM [0 The pattern of proof used to prove theorem 
[T] works equally well for some sub-elliptic convection-diffusion operators L x of Hormander 
form | Y^a=i Vi + B satisfying the following assumptions. 

• The vector fields Vi(-;y) and B(-;y) on IR n are measurable functions of (x,y) G 
W 1 x of class C°° as functions of x G M n , with all their derivatives bounded 
uniformly in y G 

• The vector fields B(-,y) satisfy hypothesis (ii) of theorem [H 

• Set A = {0} U UfcLi{l> • • • iP} k y se ^ ^[0] = an d = Vi, for i = 1 . . . ,p, and define 
inductively the vector fields VLi, for a G A, setting 

V[ ai ] = [V [a ],Vi\, for % = 1, . . . ,p. 

Given a positive integer £, write At for {0} U Ujt=i{l> • • • >P} k - We suppose there 
exists an i such that 

(2.4) inf { {V[«](x,v)>Z) 2 ; & v) e k" x m;} > 0. 

a€A e 

Indeed, theorem 2 of |Kus03| provides a uniform in a G At estimate of 

sup \V [a] P t d " b f\(x,y) 

in terms of ||/||oo) a t a fixed time t > 0, similar to (12. ip . The y-uniform equi-continuity of 
P t dlvb f follows from condition (12. 4 p as above. The remainder of the proof of theorem [T] does 
not need any change. 

This kind of improvement of theorem [T] seems difficult to get by the analytical methods 
developped in the works mentionned in the introduction. 



3. Multiple coagulations, fragmentation and scattering 

3.1. Multiple coagulations. It is clear from the proof of theorem [1] that it can be extended 
to a multiple coagulation framework. Write K2 for the above binary coagulation kernel K 
and let K n be for all n ^ 3 a symmetric non-negative valued function on • Set 

K n (n,--- ,n) = ^ J {5 yi+ ... +yn - (5 yi + ■ ■ ■ + 8 yn )}K n (y lr ■ ■ ,y n ) //(dyi) • ■ • n{dy n ) 

The spatial coagulation equation takes in that general framework the form 
(3.1) dtUt = L* x fi t + ^K n (fi t , . . . ,fi t ). 

THEOREM 6. Suppose the coagulation kernels K n are uniformly bounded above by a constant 
independent of n. Suppose hypotheses (i) and (ii) of theorem^ hold. Then equation (13. ip 
has a unique maximal solution in N ', defined for all times if \\fJ>o\\ is small enough; it is 
non-negative if [j,q is non-negative. 

PROOF - The proof of this statement is similar to the proof of theorem [U noticing that we 
have ^ C n u under condition (11.51) . and introducing, for any measure \x on M^j_, the 
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measures K£((-i, ■■■,[■!>) and K~(p,, . . . , fj) defined by the formulas 

K~{^...,p){dz) 

Denote by M an upper bound for all the coagulation kernels, and write m for i/o(M+). 
To prove that the unique maximal solution is defined for all times for a small enough 
initial condition, remark that equation (13. ip implies the inequality 

\\lh\\<e ^ ||+^(-— + --— ) / e « n»s\\ n ds, 

so \\pt\\ is bounded above by the solution of the differential equation 

it = -ezt + M(e Czt - 1 - Cz t ) + z t (e mzt - l), z = \\p \\. 
Its solution is decreasing and defined for all times provided Zq is small enough. t> 

3.2. Coagulation, fragmentation and scattering. More realistic models of chemical 
reactions have fragmentation and scattering terms. The fragmentation term accounts for 
the spontaneous breakage of particles into smaller particles; it is defined by the following 
formula, where p, here and below, is used as a generic measure on M^_, 

<8(p)(dz) := J {F(y;dz)-5 y (dz)}B(y)p(dy). 

A particle of mass y splits into smaller particles at rate B(y), the result of the splitting 
being given by the non-negative measure F(y ; dz) with support in (0, y) and total mass 
j z F(y ; dz) = y. The scattering term accounts for the fact that real particles have a mass 
no bigger than some maximum mass y$, in a given situation, so if coagulation forms a 
particle with mass bigger than yo> this particle breaks instantaneously into smaller particles. 
The scattering term (3 is determined by a scattering kernel S from (y ,2y ) to the set of 
non-negative measures on and the formula 

&{p,p)(dz) := / {S(y + y';dz) - 5 y (dz) - 5 y ,(dz)}l y+y/>yo K(y,y') p(dy)p(dy'). 

This term was first introduced in |FR00j . and studied further in |Wal02] and |Wal05] . The 
coagulation term is changed in this model into 

K(p, p) = J {S y+y/ -b y - 5 yl }l y+y ^y Q K(y, y') p{dy) p(dy'). 

The proof of the following theorem is analogue to the proof of theorem [TJ 

THEOREM 7. Let y > be an upper bound for the size of the particles in the system. 
Suppose the following conditions hold. 

(i) There exists some positive constant a, {3 such that aid ^ a(x, y) ^ (3ld, for all 
(x,y)eR n x(0,y o ). 

(ii) There exists a positive constant e such that the vector field b(-,y) has divergence 
bounded below by e, for all y G (0, yo). 



i / "2/1 H hy™ 

(dz)K n (y 1 , ■■■ ,y n ) p{dy x ) ■ ■ ■ p(dy n ), 
3— jry ( / K(yi, ■ ■ ■ ,y n -i,z) p{dyi) ■ ■ ■ p{dy n ^p(dz). 



[n — 
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(iii) The coagulation and breakage rates K, B are bounded, and the fragmentation kernel 
F(y;dz) has the form f(y;z)u (dz), with ^^(l + ||f||oo) < £■ 

(iv) We have 

(3.2) J S{y + y'; dz) v (dy)v {dy') ^ C v Q {dz) 

for some positive constant C. 
Given fio £ M , the equation 

fi t = P t fio + \ Pt- s (K(n s ,n s ) + *B(// S ) + &(iJ, s ,fj, s ))ds, te[0,T] 
Jo 

has a unique solution in M continuous with respect to t, and defined for all times provided 
ll/ioll is small enough; it is non-negative if fio is non-negative. 

Condition (iv) will for instance hold if v = LeB|( 0jJ/0 ) and the measures S(a ; dz) have a 
uniformly bounded density with respect to Lebesgue measure on (0,yo), for a £ (yo,2yo). 
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